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Abstract : 

We study complexified elliptic Calogero-Moser integrable systems. We determine the value of the 
potential at isolated extrema, as a function of the modular parameter of the torus on which the integrable 
system lives. We calculate the extrema for low rank B,C,D root systems using a mix of analytical and 
numerical tools. For so( 5) we find convincing evidence that the extrema constitute a vector valued 
modular form for the Fo(4) congruence subgroup of the modular group. For so(7) and so(8), the extrema 
split into two sets. One set contains extrema that make up vector valued modular forms for congruence 
subgroups (namely Fo(4), T(2) and r(3)), and a second set contains extrema that exhibit monodromies 
around points in the interior of the fundamental domain. The former set can be described analytically, 
while for the latter, we provide an analytic value for the point of monodromy for so(8), as well as extensive 
numerical predictions for the Fourier coefficients of the extrema. Our results on the extrema provide a 
rationale for integrality properties observed in integrable models, and embed these into the theory of 
vector valued modular forms. Moreover, using the data we gather on the modularity of complexified 
integrable system extrema, we analyse the massive vacua of mass deformed A/” = 4 supersymmetric 
Yang-Mills theories with low rank gauge group of type B , C and D. We map out their transformation 
properties under the infrared electric-magnetic duality group as well as under triality for J\f = 1* with 
gauge algebra so(8). We compare the exact massive vacua on IR 3 x S 1 to those found in a semi-classical 
analysis. We identify several intriguing features of the quantum gauge theories. 
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1 Introduction 

Four-dinrensional gauge theories accurately describe forces of nature. Since solving them is hard, we may 
revert to studying supersymmetric four-dinrensional gauge theories, in which the power of holomorphy 
lends a helping hand. Twenty years ago, we realised how to solve for the low-energy effective action on 
the Coulomb branch of Af = 2 gauge theories in four dimensions m- The solution techniques were soon 
recognised to lie close to those studied in integrable systems [3 SI- It is the bridge between integrable 
models and supersymmetric gauge theories that we will further explore in this paper. We also attempt 
to reinforce both sides separately, and present results in a manner such that the contributions to these 
two domains may be read independently. 

The link between supersymmetric gauge theories and integrable systems was useful in writing down 
the low-energy effective action for AT = 2 * gauge theory, namely Af = 4 super Yang-Mills theory with 
gauge group G, broken to Af = 2 supersymmetry by adding a mass term for one hypermultiplet. For 
the gauge group G = SU(N) this program was completed in terms of a Hitchin integrable system with 
SL(N, C) bundle over a torus with puncture [F- The associated elliptic Calogero-Moser system permits 
generalisations to any root system, and allows for twists, which were used to provide Seiberg-Witten 
curves and differentials for Af = 2* theory with general gauge group G [ 6 ]. The generalisation was non¬ 
trivial since the elegant technique of lifting to M-theory [7] is difficult to implement in the presence of 
orientifold planes (see e.g. M), while the relevant generalised Hitchin integrable system has a gauge 
group which is related to the gauge group of the Yang-Mills theory in an intricate manner m- For a 
review of part of the history, see the lectures P33- 

We will be interested in breaking supersymmetry further, from Af = 2 to Af = 1 by adding another 
mass term for the remaining chiral multiplet (providing us with three massive chiral multiplets of arbitrary 
mass). We will study this Af = 1* gauge theory with generic gauge group G. With Af = 1 supersym¬ 
metry, we hope to calculate the effective superpotential W at low energies exactly. For an adjoint mass 
deformation from Af = 2 to Af = 1 this was done in the original work [I] in certain cases. For Af = 1* 
and gauge group G = SU(N), the exact superpotential was proposed in [12] following the techniques of 


2 



















mm- The superpotential is the potential of the complexified elliptic Calogero-Moser integrable system 
associated to the root lattice of type A N _ i. In m the exact superpotential for AT = 1 * with more 
general gauge algebra was argued to be the potential of the twisted elliptic Calogero-Moser system with 
root lattice associated to the Lie algebra of the gauge group G. See [TS] for further generalizations to 
Af = 1* theories with twisted boundary conditions on R 2 3 x S' 1 . 

In this paper, we wish to analyse the proposed exact superpotential in more detail. This involves a 
study of the properties of the isolated extrema of the complexified and twisted elliptic Calogero-Moser 
integrable system. The results are of independent interest, and we have therefore dedicated a first part 
of this paper to the study of the integrable systems per se. 

The paper is structured as follows. In section[2j we review the relevant elliptic Calogero-Moser models. 
We pause to demonstrate a Langlands duality between the B and C type integrable systems. We then 
analyse the isolated extrema of the complexified potential of low rank integrable systems of B 1 C and D 
type, and their modular properties. We observe the strong connection to vector valued modular forms. 
The latter in turn provide a natural backdrop for integrality properties of integrable systems (see e.g. 
unnnuiH]). Section [2] is the technical heart of the paper, and we will lay bare many properties of the 
vector valued modular forms, using a combination of analytical work and extensive numerics. We will 
analytically describe the potential in certain classes of extrema. We also find sets of extrema that exhibit 
a monodromy in the interior of the fundamental domain. In these cases we are able to calculate the 
monodromy, as well as to provide extensive numerical data for the integer valued coefficients describing 
the value of the potential at the extrema. 

Finally, in section [3j we reinterpret the results we obtained in terms of the physics of massive vacua of 
AT = 1* theories. We compare our results for the quantum theory on R 3 x S ' 1 to semi-classical results for 
massive vacua and discuss electric-magnetic duality properties in the infrared under the modular group 
as well as the Hecke group. For so( 8 ), we also detail the action of the global triality symmetry on the 
massive vacua. We will encounter several interesting phenomena. We conclude in section [4] and argue 
that we have only scratched the surface of a broad field of open problems. 

2 Elliptic Integrable Systems and Modularity 

It is interesting to identify and study dynamical systems that are integrable. Often they form solvable 
subsectors of more complicated theories of even more physical interest. There exist one-dimensional 
models of particles with interactions that are integrable, and the Calogero-Moser models of our interest 
are one such class nm [ini up. These models are associated to root systems of Lie algebras (amongst 
others). See e.g. [22], [53] for a review. Integrable systems are also known to have certain integrality 
properties. Namely, their minimal energy, frequencies of small oscillations as well as eigenvalues of Lax 
matrices are often expressible in terms of a series of integers mi na mi- 

In this section, we study properties of (twisted) elliptic Calogero-Moser systems. We analyse the 
complexified model, defined on a torus with modular parameter r. In particular, we examine the extrema 
of the complexified potential, and exhibit their curious characteristics. 

2.1 The Elliptic Calogero-Moser models 

The member of the pyramid of Calogero-Moser integrable systems we concentrate on is the elliptic 
Calogero-Moser model. We concentrate on the models associated with a root system A, as well as their 
twisted counterparts. These models have a Hamiltonian with rank r variables, with canonical kinetic 
term, and a potential of the form: 

Va = g^2p{a(X)-u} i,w 2 ), (2.1) 

a£ A 

where p is the Weierstrass elliptic function on a torus with periods 2uq and 2 w 2 and g is a coupling 
constant. We choose the half-periods such that the imaginary part of the modular parameter r = W 2 /W 1 
is positive)^] The vector X lives in the space dual to the root lattice of rank r and the sum in the potential 
is over all the roots a of the root system AQ The model is integrable for all Lie algebra root systems. 
The twisted elliptic Calogero-Moser model is defined in terms of twisted Weierstrass functions: 

k 

Pn(x;u 1 ,w 2 ) = V' p{x + -2wi;wi,w 2 ) , (2.2) 

' n 

which are summed over shifts by fractions of periods (thus in effect modifying that period). We have a 
twisted elliptic Calogero-Moser model for all non-simply laced root systems and the value of n is then 

2 See appendix IbI for more on our conventions for elliptic functions. 

3 See appendix|7\ for our conventions and a compendium of properties of Lie algebras and Lie groups. 
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given by the ratio of the length squared of the long versus the short roots. We will be interested in the 
twisted elliptic Calogero-Moser model with potential: 

Va ,tw = 91^2 p(ai(X)-,wi,u> 2 ) +g s ^2 Pn(«sW;wi,W 2 ), (2.3) 

aieA; o.eA, 

where cq denote the long and a s the short roots in the root system A = A; U A s , and gi and g s are two 
coupling constants. We will concentrate on the root systems A r , B ri C r and D r corresponding to the 
classical algebras su(r + 1), so(2r + 1), sp{2r) and so(2r). We allow complex values for the components 
of the vector A' (i.e. X £ C r ). 

The symmetries of the potential 

Let us discuss in detail the symmetries of the twisted elliptic Calogero model that act on the set of variables 
X. We first observe that the Weyl group action leaves invariant the scalar product a(X) = (a, X) and that 
the root system is Weyl invariant [(] This implies that the Weyl group action on X leaves the potential 
invariant. Secondly, we note that the outer automorphisms of the Lie algebra, which correspond to 
symmetries of the Dynkin diagram, also leave the set of roots and the scalar product invariant. Therefore, 
outer automorphisms as well form a symmetry of the model. 

Moreover, the periodicities of the model in the two directions of the torus are as follows. By the 
definition of the dual weight, or co-weight lattice, we have that a(A v ) £ Z for all roots a. This implies 
that shifts of X by 2oj 2 P v , namely shifts by periods times co-weights, leave the potential invariant. 

To discuss the periodicity in the w\ direction, we concentrate for simplicity on the algebras A 7 B , C 
and D , and normalize their long roots to have length squared two. We then have that for a long root 
cq and a weight A, the equation (cq, A) £ Z holds while for a short root a s of the B or C algebras we 
have (a s ,A) £ |Z, for all weights A. As a consequence, the periodicity in the (twisted) ui\ direction is 
the lattice 2uq P where P is the weight lattice. The group of all symmetries is a semi-direct product of 
the lattice shifts, the Weyl group as well as the outer automorphism group. 


2.2 Langlands Duality 


Beyond the many features of these integrable systems already discussed in the literature, the first sup¬ 
plementary property that will be pertinent to our study of isolated extrema, is their behaviour under an 
inversion of the modular parameter r. We therefore briefly digress in this subsection to discuss a few 
of the details of the duality. Models associated to simply laced Lie algebras map to themselves under 
the modular S-transformation S : r —> —Xjr. This is easily confirmed using the transformation rule 


(B.3) of the Weierstrass p function under modular transformations. We do have a non-trivial Langlands 


or short-long root duality between the twisted elliptic Calogero-Moser model of B-type and the twisted 
model of C-type. In order to exhibit the duality, we make the potential for the (twisted) B r = so(2r + 1) 
theory more explicit]^] 


Vb = bi 

y p(Xi - Xj-ljJ!,^) + p(Xi + Xj- 7 U 1 ,W 2 ) 

+ b s 

r 

y p{x£ UJ 1 , U} 2 ) + p(xi + uq; uq, ui 2 ) 


i<j 


_i =1 


and for the C r = sp(2r) theory as well: 


Vc = 


i<j 


p{yi - yj',u[,uj 2 ) + p(yi + yf,u[,u' 2 ) + p(yi - y 3 + (,) + p(y t + y 3 + a/,; a/,, w(,) 


+ci'22p(2y i -,u}'[,u/ 2 ). 


(2.4) 


We have chosen a standard parameterisation of the vector X as well as the root systems, and we have 
assigned half-periods oji to the B-system and w' to the C-system. We have also made explicit the twisted 
Weierstrass functions p 2 with twisting index 2, which is the ratio of lengths squared of the long and short 


roots. To demonstrate the duality between these models, we use the elliptic function identities (B.5| to 
manipulate the so(2r + 1) potential such that it becomes of the form of the sp(2r) potential: 


V B 


bi 


Y, p{xi - x 3 ; 2w 2 , -wi) + p(xi - Xj + 2 w 2 ; 2w 2 , -uq) + p(aq + Xf, 2u 2 , -uq)+ 

i<j 


4 We mostly follow m for our conventions on Lie algebras. See also appendix [a] for the definitions of the different lattices 
discussed hereafter. 

5 For the non-simply laced cases, we will always work with the twisted model, and we will drop the corresponding subscript 
on the potential from now on. 
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p(xi + Xj + 2 oj 2 ; 2 w 2 , - wi ) 


+ b s y p(xi;uj 2 , —uii/2) 


7 r 2 r(r — 1) 
24w| 




Wl 


2E 2 l -- 1 - Es I - 

w 2 


_Wl_ 

2 w 2 


7 r 2 r 

6 ^ 


>w 2 


Wi 


2 £, 2 - - 


w 2 


Wi 


= bt 


y, p(xi - Xj; 2w 2 , -Wi) + p(xi - Xj + 2w 2 ; 2w 2 , -wi) + p(x* + Xjj 2w 2 , -Wi)+ 


i<j 


p(xi + Xj + 2oj 2 ; 2w 2 , — Wi) 


+ I2^( 2r6s+ r ( r _ i) 6 ') 


+ 6 s ^p(xj;o; 2 ,-Wi/ 2 ) 


2£ 2 ( 2 —-£ 2 ( — 
Ui ) \UJi 


(2.5) 


In the last equality, we used the modular transformation rule (C.3) for a combination of second Eisenstein 
series. We observe that the end result (2.5) can be identified with the C r potential (2.4), provided we 
match parameters as follows: 


= 2w 2 u}' 2 = —cj i yi = Xi 

c s =bi ci = 46 s , 


( 2 . 6 ) 


and we allow for a r-dependent shift of the potential that invokes the second Eisenstein series E 2 . These 
identifications imply a duality (which we will denote S 2 ) between the modular parameters of the B and 
C-type integrable systems: 


r 


B _ 


1 

2 T C 


(2.7) 


In the following, we will be interested in B and C models in which the ratio of the long to short root 
coupling constants is equal to two, i.e. we put bi = b = 2 b s and c; = c = 2c s j^] It is important that 
this relation is compatible with the duality map (2.6). We rewrite the identity of the potentials for this 
specific ratio of parameters: 


y p(xi - Xj;u)\,L0 2 ) + p(xi + Xj;wi,u} 2 ) + \ ( y p(xi;wi,w 2 ) + p(x* + wi;wi,w 2 ) j = 

i<j \ i / 


y p{xi - Xj; 2 uj 2 , -ui i) + p(xi - Xj + 2u> 2 ; 2 w 2 , -uq) + p(x* + Xj; 2w 2 , -uq) + 

i<j 

p(xi + Xj + 2 w 2 ; 2w 2 , —uq) + 2 ^ ) p[2xt; 2w 2 , —uq) + 


12 wf 


,w 2 


Wi 


2 E 2 2 — - — 


w 2 


Wi 


and the integrable system duality can be summarised as: 


1 


Xi 1 


V B (xi,T) = — 2 Vc I ) + — [2E 2 (2r) - E 2 (r)\ , 


7 r 2 r 2 


( 2 . 8 ) 


when we use the rescaling (B.2). The duality may be viewed as a standard Langlands duality. We went 
through its detailed derivation since the r-dependent shift in the duality transformation (2.8) is important 
for later purposes. 


2.2.1 Langlands duality at rank two 

There is a further special case of low rank which is of particular interest to us in the following. The B 
and C type Lie algebras of rank two are identical: so(5) = sp{ 4). If we apply the duality of B and C type 
potentials to this special case, we derive that the following transformations leave the potential invariant: 

= 2u 2 uj 2 = —oq c' = 2b 

x 2 — x\ = 2xi x\ + x' 2 = 2x 2 . (2-9) 

6 Various particular choices of parameters and observables that we make in section [3] are motivated by the gauge theory 
applications that we will discuss in section |3l It is also of interest to study the integrable systems more generally. 
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If we parameterise the potential in terms of the modular parameter r = w 2 /wi, the duality transformation 
for so(5) reads: 


Ko(5) (*^1 > 7~) 2^-2 ^°(5) ^ 


X\ + X 2 X\ — X 2 1 
2 T ’ 2t ’~2T 


47 T J 


[2E 2 (2t) - P 2 (r)] . 


( 2 . 10 ) 


In summary, we derived a Langlands duality between B and C type (twisted) elliptic Calogero-Moser 


models. The resulting identities captured in equations (2.8) and (2.10) and the shifts appearing in these 


duality transformations will be useful. We return to the more general discussion of the integrable systems, 
and in particular their extrema. 


2.3 Integrable Models at Extrema 

There have been many studies of classical integrable models at equilibrium. These have uncovered 
remarkable properties, like the integrality of the minimum of the potential and of the frequencies of small 
oscillations around the minimum, amongst others (see e.g. [TOUT?! [T8] ). We will analyse the potential of 
certain elliptic integrable systems evaluated at generalised equilibrium positions. We show that they give 
rise to interesting vector valued modular forms as well as more general non-analytic modular vectors. 
Modularity provides a more conceptual way of understanding the integrality properties of the integrable 
system. This rationale then continues to hold for the integrable systems that can be obtained from the 
elliptic Caloger-Moser systems by limiting procedures (e.g. the trigonometric models). Thus, studying 
elliptic integrable systems, depending on a modular parameter, is found to have an additional pay-off. 

It is known that A-type integrable systems often have simpler properties than do the integrable 
systems associated with other root systems. As a relevant example, let us quote the fact that the 
(real) Calogero-Moser (Sutherland) system with trigonometric potential of A-type has equally spaced 
equilibrium positions along the real axis, while the B , C, D-type potentials have minima associated to 
zeroes of Jacobi polynomials in, which satisfy known relations [55], but are not known explicitly in 
general. The elliptic Calogero-Moser systems that we examine show a similar dichotomy. Extrema of the 
(complex) elliptic A-model are equally spaced. This fact leads to relatively easily constructable values 
for the potential at extrema, for any rank umm- For the B 1 C, D -type models that we study in 
this paper, much less is known, and we need to combine numerical searches with analytic approaches to 
determine the extremal values of the potential, for low rank cases. 

To be more precise, we will be interested in extrema of the complexified potential, satisfying^] 


d Xi v(x j ) = o w, 


and we moreover demand that at the extremum (2.11) the function 


5 > X W(W )| 2 

2=1 


( 2 . 11 ) 


( 2 . 12 ) 


not posses any flat directions]^] 

Recall that the group of symmetries acting on the variables X were a lattice group of translations, 
the Weyl group as well as the outer automorphisms of the Lie algebra. Using these symmetries, we will 
introduce a notion of equivalence on the variables X. We will consider the vector X to be identified by 
the periodicities of the model. The periodicity in the u> ± direction is given by the weight lattice P, while 
in the w 2 direction it is the co-weight lattice P v . Furthermore, we will consider extrema that are related 
by the action of the Weyl group of the Lie algebra to be equivalent. By contrast, outer automorphisms 
are taken to be global symmetries of the problem. When the global symmetry group is broken by a given 
extremum, the global symmetries will generate a set of degenerate extrema. 


2.4 The Case A r = su(r + 1) 

The extrema of the elliptic Calogero-Moser model of type A r have been studied in great detail, mostly 
in the context of supersymmetric gauge theory dynamics (see e.g. ill [H] [26]). Firstly, we remark that in 
this case, the equivalence relations that follow from the periodicity of the potential as well as the Weyl 
symmetry group of the Lie algebra are straightforwardly implemented. We use the parameterisation 
of simple roots in terms of orthogonal vectors at = ej — Ci+i, and the fundamental weights then read 
7 Tj = Xq=i e B with weight lattice spanned by the vectors e,. We can parameterise the coordinates of our 

7 This will correspond, in section [ 3 ] to a supersymmetric vacuum in the J\f = 1* gauge theory, where the effective 
superpotential W is identified with the potential V of the integrable system. 

8 This condition implies that the vacuum is massive in the supersymmetric gauge theory. We briefly comment on massless 
vacua later on. 
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integrable system by a vector Xjej living in the dual to the root space (and = 5B). The Weyl 

group S n acts by permuting the components Xj. We can shift one of the components Xj to zero by 
convention. The equivalence under shifts by fundamental weights is identical to the toroidal periodicity 
relations for the individual coordinates Xj. The inequivalent extrema of the su(n) potential (satisfying 
the additional condition ( 2 . 12 ) of non-flatness) are then argued to correspond one-to-one to sublattices 
of order n of the torus with modular parameter r mm- These extrema are classified by two integers 
p and k satisfying that p is a divisor of n and k £ {0,1, — 1}. The number of extrema is equal to 
the sum of the divisors of n. The Z 2 outer automorphism of A r>1 acts trivially on the minima, since it 
acts by permutation, combined with a sign flip for all Xj , which leaves a sublattice ankered at the origin 
invariant. 

The value of the potential at one of these extrema is (with a given choice of coupling constant): 


Va^At) = 


n 3 („ . . P „ (P k 
-(E 2 (r)--E 2 (-r+- 


(2.13) 


Under the <SX(2,Z) action on the torus modular parameter r, the sublattices of order n of the torus are 
permuted into each other (in a way that depends intricately on the integer n). The permutation of the 


sublattices also entails the permutation of the values (2.13) at these extrema under SL{2,T,). The list of 
extremal values of the elliptic Calogero-Moser model therefore form a vector valued modular form (see 
e.g. [menEss ) of weight two under the group <SX(2,Z). The associated representation of the modular 
group is a representation in terms of permutations specified by the SL( 2, Z) action on sublattices of order 
n. One can identify a subgroup of the modular group under which a given component of the vector-valued 
modular form is invariant, and then use minimal data to fix it l30j . 

In summary, the extrema of the Calogero-Moser model of type A r = su(r + 1) are under analytic 
control. The positioning of the extrema can be expressed linearly in terms of the periods of the model, 
and the vector valued modular form of extremal values for the potential has an automorphy factor that 
can be characterised by sublattice permutation properties. The extremal values are generalised Eisenstein 
series of weight two under congruence subgroups of the modular group. 


2.5 The B , C , D Models 

For other algebras, we are at the moment only able to study low rank cases. From the analysis, it is 
clear that crucial simplifying properties of the A r case are absent. Nevertheless, generic features of the 
A r case persist in a subclass of extrema, in that we find vector-valued modular forms as extremal values 
for the potential. We also find a class of extremal values that exhibit new features. 

To describe in detail which extrema are considered to be equivalent, we must discuss the equivalence 
relations that we mod out by for the B , C and D root systems individually. 

D r = so(2r ) 

For the D r case, we can parameterise the roots as — e^+i (for i £ {1,2,... ,r — 1}) and a r = 

e r -i + e r . We put X = Xje? and imply that the relation e^e- 7 ) = 6, J holds. The equivalence of the 
vector X under shifts proportional to the weight lattice implies that each variable Xj lives on a torus 
with modular parameter r. It moreover identifies the vector X with the vector X shifted by a half-period 
in each variable simultaneously. The Weyl group is W(so(2r)) = S r k Z 2 ~ , and acts by permutation 
of the components Xj, as well as the sign change of an even number of them. The outer automorphism 
group (for r/ 4) is equal to Z 2 and acts as X r —► —X r . For r = 4, the global symmetry group is S 3 
triality. 

B r = so(2r + 1) 

For B r , the roots are on = e* — ej+i (for i £ {1, 2,..., r — 1}) and a r = e r . We recall that the periodicity 
is the weight lattice in the uq direction (due to the twist), and the co-weight lattice in the w 2 direction. 
Thus, we can shift components of the vector X = Xje 3 by periods, or all components simultaneously by 
a half period in the u> 1 direction. In the w 2 direction, we allow shifts of the individual components by 
periods. The Weyl group acts by combinations of permutations and any sign flip of the coordinates. 

C r = sp{2r) 

The roots are a, = (e, — e i+ \)/\f2 (for i £ {1, 2,..., r — 1}) and a r = 1 / 26 ^ We can shift components Xj 
of X = y/2Xje^ by half-periods in the uq direction, while in the w 2 direction, we can allow shifts by any 
period, as well as a half-period shift of all Xj simultaneously. The Weyl group allows any permutation 

lJ By our conventions, we normalise the long roots such that they have length squared two. 
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and sign flip of the coordinates. The equivalence relations and symmetries in the B , C and D cases, 
beyond permutation symmetries and toroidal periodicity, are summarized in the table: 


B r 

Individual Xj — > — Xj 

Collective Xi — ► V + wi 

C r 

Individual Xj — > — Xj and Xj — > Xj + uq 

Collective Xj — y Xj + co 2 

D r 

Even number of sign flips Xj — > — Xj 

Collective Xj — > Xj + u>i and Xj —► Xj + oj 2 

Global symmetries : Z 2 generically and S 3 for D 4 . 


Armed with this detailed knowledge about the equivalence of configurations, we programmed a nu¬ 
merical search for isolated extrema. In the following subsections, we list the results we found by root 
system. For simply laced root systems we studied the elliptic Caloger-Moser model, while results for non- 
simply laced root systems correspond to the twisted elliptic Calogero-Moser model with a coefficient for 
the short root term which is equal to one half the coefficient in front of the long root terms (as described 
below equation (2.7)). 


2.6 The Case C 2 = sp( 4) = so(5) and Vector Valued Modular Forms 

Since the root system C 2 is the first example of our series, we provide a detailed discussion. We discuss 
the positions of the isolated extrema, the series expansions relevant to the potential at these extrema, 
the action of the duality group, as well as the identification of the relevant vector valued modular forms. 


2.6.1 The positions of the extrema 

For the Lie algebra so(5) = sp( 4) we found 7 isolated extrema of the potential. We provide their 
positioning at t = i in figure [l] We have drawn in bold the positions of the extrema as well as their 
opposites, in a fundamental cell of the torusp 7 ] 

These numerical results were found using a Mathematica program, which was written around the 
built-in function FindMinimum. Careful programming augments the precision of the algorithm to at least 
two hundred digits. The most costly part of the algorithm is the random search for extrema. Indeed, 
the intricate landscape drawn by the potential can hide extrema. We gave a drawing of the position 
of the numbered extrema on the torus with modular parameter t = i. The positions of the extrema 
for other values of the modular parameter can be reached by interpolation. We have analytic control 
over a few extra properties of the extrema. E.g. if we follow extremum 1 to r = ioo, we find that 
the equilibrium positions are given by ^ arccos(±l/-\/3) where ±l/v^3 are the zeroes of the Jacobi 
polynomial P^ 0,0 ^. The first extremum, which we label 1, lies on the real axis and is the equilibrium 
position of the real integrable system. The extremum 2 lies on the imaginary axis, while extrema 3 and 
4 are then approximately obtained by applying the transformation r —> r + 1. The extrema 5 and 6 are 
S 2 Langlands duals of extrema 3 and 4. It is easy to deduce from the potential that the positions of the 
extrema generically behave non-linearly as a function of r. 


2.6.2 Series expansions of the extrema 

By numerically evaluating the extrema of the potential for a range of values of the modular parameter r, 
we are able to write the extrema as an expansion in terms of a power of the modular parameter q = e 2mT . 
The extremal values can be written in terms of the series: 

A 0 {q) = ^+q + q 2 +4q 3 + q 4 + 6q 5 +4q 6 + 8q 7 + q S + 13q 9 + 6q W + 12Q 11 

+4q 12 + Uq 13 + 8q 14 + ... (2.14) 

A x {q) = 1 + 48g + 828 q 2 + 8064q 3 + 109890g 4 + 1451520g 5 + 11198088<? 6 + 141212160<j 7 
+1666682811g s + 9413050176g 9 + 145022264892 9 10 + 1838450006784< ? 11 
+11103941590326g 12 + 138638111404032g 13 + ... 

A 2 {q) = 2 + 48 q + 576 q 2 + 9792 q 3 + 99576 q 4 + 743904g 5 + 13146624g 6 + 115737984g 7 

10 We have indicated reflections over other half-periods in grey, to illustrate that the minima are close to forming sublattice 
structures. 
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Figure 1: Extrema at r = * for the Lie algebra so(5) 



Extremum 7 



+1015727364g 8 + 14338442448g 9 + 102050482176g 10 + 935515738944g n 
+12532363069968g 12 + 122390111091744g 13 + ... 

A 3 (q) = ^ + Iq + 541g 2 + 24508g 3 + 939669g 4 + 19944842g 5 + 764752180g 6 

+21016537080g 7 + 905672825157g® + 38827071780859g 9 + 827503353279726g 10 + ... 
Ai(q) = 1 + 148g + 7446g 2 + 154344g 3 + 5100349g 4 + 352720380g 5 + 10627587582g 6 
+166124184888g 7 + 5419843397586g® + 294399334337124g 9 + ... 

A 5 (q) = - — + 29g + 431g 2 + 80468g 3 - 231081g 4 + 94846414g 5 + 1301490428g 6 

+90560563752g 7 - 529100109849g® + 93349951292249g s + ... . (2.15) 


The integer coefficients have been determined up to an accuracy of at least 10~ 6 . For the first order 
terms, the accuracy can be up to 10 -2no . In terms of these series, the potential in extremum number 1, 
on the real axis is (with a given choice of normalisation): 


Vi = 144^3 (J^) ■ 


The potential in the other extrema are: 


and 


(2.16) 


V 2 = -12tt 2 Q,4o(g) + (2q) 1 '' 3 A 1 (q/9) + (2q) 2/3 A 2 (q/9)^j 

v 3 = -127T 2 ^Ao(q) + (2 q ) 1 / 3 e 2m / 3 A 1 (q/9) + (2q) 2 / 3 e 4 ^ 3 A 2 (q/9)^) 

V 4 = - 12 ^ 2 ^ 0 (g) + (2g) 1 / 3 e 4 "/ 3 ^ 1 (g/9) + (2g) 2 / 3 e 27 ri / 3 A 2 (g/9)^ , 


(2.17) 


^ 5,6 = 


V 7 = 



(2.18) 


The growth properties of these series, as well as the fact that we are dealing with a physical system living 
on a torus suggests turning these numerical data into an analytic understanding, based on the theory of 
modular forms. In the following, we show that this is possible for the rank 2 root system B- 2 - 


2.6.3 The extrema as modular forms of the Hecke group and the T 0 (4) subgroup 


We need to introduce a few groups related to the modular group. We already noted the duality transform 
for the B,C -type twisted Calogero-Moser system under the map S 2 : r —>• — 1/(2 t) (see equation (2.7)). 
For the so( 5) Lie algebra, which is identical to the sp( 4) Lie algebra, this transformation maps the 
integrable system to itself (up to a r dependent shift of the potential and an overall factor - see equation 
(2.10)). The map T : t —» r + 1 also maps the integrable system to itself. Together, these transformations 
generate the action of a Hecke group dubbed T*(2) on the modular parameter r. This group contains a 
subgroup To(4) which is a congruence subgroup of the modular group SL( 2,Z). Generators of the group 
Tq (4) can be chosen to be the 2x2 matrices: 


T 

U 



(2.19) 


The action of these matrices on r coincides with the action of the elements T and U = S 3 T 2 S 2 of the 
Hecke group. For more information on Hecke groups and associated modular forms see e.g. the lectures 


The extremal values of the potential may therefore form a vector valued modular form with respect 
to the Hecke group T*(2), and as a consequence also with respect to the congruence subgroup r 0 (4) of 
the modular group SL( 2,Z), since we expect extrema to be at most permuted and/or rescaled under the 
group. Here, we assume analyticity in the interior of the fundamental domain. We will mostly exploit 
the group Fq( 4) in the following, since the literature on the subject of modular forms with respect to 
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Figure 2: The diagram of the action of dualities on the extrema for £2 = so( 5). In red, we draw the 
action of Langlands ^-duality, and in green, T-duality (when the action is non-trivial). 


congruence subgroups is abundant. For starters, we determine the action of the operations T and S 2 on 
the vector V) of extremal values of the twisted Calogero-Moser potential: 


T 


S 2 


/ 10 0 0 
0 0 0 1 
0 10 0 
0 0 10 
0 0 0 0 
0 0 0 0 
\ 0 0 0 0 

/ 0 1 0 0 
10 0 0 
0 0 0 0 
0 0 0 0 
0 0 10 
0 0 0 1 
\o 0 0 0 


0 0 0 \ 

0 0 0 

0 0 0 

0 0 0 

0 1 0 

1 0 0 

0 0 1 / 

0 0 -2 \ 

0 0-2 
1 0 -2 

0 1 -2 

0 0-2 
0 0-2 
0 0 - 1 / 


( 2 . 20 ) 


See figure [2] for a summary of the action of the duality group. To this information, we add the last 
column in the matrix S 2 , which originates in the shift of the potential under Langlands duality. From 
these data, we easily calculate the action of the generator U = S 2 T~ 2 S 2 on the vector valued modular 
form: 


U 


/ 0 0 
0 1 
0 0 
0 0 
0 0 
1 0 
\ 0 0 


0 0 1 
0 0 0 
1 0 0 
0 1 0 
0 0 0 
0 0 0 
0 0 0 


0 0 \ 
0 0 
0 0 
0 0 
1 0 
0 0 
0 1 j 


( 2 . 21 ) 


We thus find the action of r 0 (4) on the vector valued modular form, and we observe the following pattern: 
there is one entry (the seventh) which is an ordinary modular form of weight 2 under r 0 (4), and there 
are two sets of three components (namely {2,3,4} and {1,5,6}) that mix under r 0 (4). Thus, our vector 
valued modular form of dimension seven splits into a singlet and a sextuplet. Concentrating on the 
ordinary modular form of weight 2, we have that it is a linear combination of Eisenstein series £ 2 ,at 
defined by: 


E 2 ,n{t) = E 2 (t) — NE 2 (Nt) . 


( 2 . 22 ) 


Indeed, the dimension of the space A42(r 0 (4)) of modular forms of T 0 (4) is two, and it is spanned by 
£ 2,2 and £ 2 , 4 - We thus only need two Fourier coefficients to fix the entire modular form, and we find 
that: 


Mq) = -^E 2 ,2(T) = ^(0j + 6j)(T) (2.23) 

Vr = y(^ + £l)(r). (2.24) 


We then have a slew of consistency checks on all the other integers that we determined numerically (see 
(2.14)). These thirteen checks work out. We do therefore claim that the result (2.24) is exact. This is a 
simple example illustrating our methodology. 

Next, we consider the triplet consisting of the components {2,3,4}. We find three eigenvectors of 
T, with eigenvalues corresponding to the cubic roots of unity. The eigenvector with eigenvalue 1 is also 
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mapped to itself under the U transformation, and forms again a modular form of weight 2 under To (4). 
It is indeed proportional to E 22 : 

V 2 +V 3 + V 4 = -2Tr 2 (0j+0j)(T). (2.25) 

The other two eigenvectors, we raise to the power three, such that they become invariant under the T- 
transformation. These forms belong to the space Al6(r 0 (4)) of weight six modular forms. The dimension 
of this vector space is 4 (see theorem 3.5.1 in [32; with g = £2 = £3 = 0 and = 3), and it consists of 
three Eisenstein series, and one cusp form. A basis for these vector spaces is given by: 


E\ 


(2.26) 

El 

= -2H2 £ «< 2r > 

(2.27) 

El 

= -^ a (4r) 

(2.28) 

S e 

= o(« 2 ) 12 , 

(2.29) 


where Eq is the Eisenstein series of weight six, and 77 is the ^-function, also recorded in appendix [C| 
We need four coefficients to fix the eigenvectors in terms of this basis and we find (using the notation 
ui 3 = exp(27ri/3)) : 

{V 2 + W 3 E 3 + w|C 4 ) 3 = — 233287 r 6 (Eg - E% - 2 S 6 ) (2.30) 

(V 2 + u>lV 3 + uj 3 V 4 ) 3 = - 23328 tt 6 (£ , 6 1 -E 2 + 2S 6 ). (2.31) 

The consistency checks using the numerics work out. 

For the second triplet, we diagonalise U first, and proceed very analogously as above, except that we 
have to take a higher power for the second combination to find a modular form of weight 12 with respect 
to To(4). We find the relations: 

(V) + W3V5 + uj 3 Vfj) 3 + (V) + u > 3 V5 + uj 3 V §) 3 = 58 327r 6 (£’g (g) — 64Eg) 

((Vi + w 3 E 5 + w 3 V 6 ) 3 - (Ei + W 3 V 5 + oj 3 V 6 ) 3 ) 2 = 1360488967 r 12 77 (g ) 24 . (2.32) 

Note that the sum of all potentials is necessarily a modular form with weight 2 of To(4). Indeed, this 
sum is equal to 1127r 2 A 0 (g) (as follows from the identity A 5 (q/27) + A 3 (q/27) = | A 0 (q)). 

2.6.4 A remark on a manifold of extrema 

There are also branches of extrema, namely, non-isolated extrema. These too, we expect to behave well 
under a modular subgroup. Although this was not the focus of our investigation, we did find numerical 
evidence for a manifold of extrema at which the potential takes the To (4) covariant value —^-E 2t2 . 

Summary 

In summary, we have full analytic control over the value of the potential for all isolated extrema of the 
so(5) twisted Calogero-Moser integrable system. We have found a vector valued modular form of weight 
two of To (4), and we were able to explicitly express its seven components in terms of ordinary modular 
forms of T 0 (4). The vector valued septuplet splits into a singlet modular form and a sextuplet vector 
valued modular form. The plot will thicken at higher rank. 


2.7 The Case D 4 = so(8) and the Point of Monodromy 

At this stage, we choose to present our results on the rank four D 4 = so( 8 ) model first, since they are 
simpler than those on the non-trivial rank three cases to be presented in subsection 2.8 The so( 8 ) model 
is simply laced and we therefore expect the ordinary modular group SL( 2, Z) to play the leading role. 
The integrable system exhibits a global symmetry group S 3 that permutes the three satellite simple roots 
of the Dynkin diagram of so( 8 ). We will refer to the S 3 permutation group as triality. We turn to the 
enumeration and classification of the extrema of the potential. We found 34 extrema. These are listed 
and labelled in appendix D.l If we mod out by the global symmetry group, we are left with 20 extrema. 


The latter fall into multiplets of the duality group of size 1,3,4 and 12. We discuss these multiplets in 
the following paragraphs. 


2.7.1 The singlet 

There is a singlet under S and T duality as well as triality. It has zero potential: V) = 0. 
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2.7.2 The triplet 

There is also a triplet under the duality group, labelled {2, 3,4}, and the dualities act as: 


T = 



S = 



The relations S 2 = 1 and {ST) 3 = 1 are satisfied. We note that in these extrema, the positions belong 
to the lattice generated by wi/2 and uj 2 / 2. For this multiplet,T-duality acts geometrically. 

We would like to deduce again from the S and T matrices and from the known first coefficients of 


the series expansions (see appendix D.l) the exact expressions of the potentials in these extrema. The 
functions are expected to transform well under some congruence subgroup of the modular group. Note 
that the sum of the three functions must be a full-fledged modular form indeed, the sum V 2 (q) + 
V 3 (q) + 14 (?) vanishes. A brute force strategy leading to the identification of the appropriate congruence 
subgroup is the following. We decompose the generators of congruence subgroups in terms of a product 
of S and T operations. We evaluate the product using the representation at hand (here 3x3 matrices) 
and check whether it is trivial for every generator. 

It turns out that the subgroup T(2) acts trivially on the extremal potentials. Hence all the potentials 
V 2 , V 3 and V 4 belong to A4 2 (r(2)). This space has dimension 2, and it is the set of linear combinations of 
the three Eisenstein functions associated to the three vectors of order 2 in (Z 2 ) 2 which have the property 
that the sum of the three coefficients vanishes. (See appendix [C] for details and conventions). Matching 
a few coefficients, we find that 


V 2 

^3 

v 4 


12 ^2G 2j2 
12 (— G 2j2 
12 ( —G 2j2 


-G 2i2 
— g 2j2 

+ 2G 2 9 


~G 2 ,2 

+ 2 G 2j2 
— G 2j2 


This can also be written in terms of the Weierstrass p function : 

um = 3 ( 2 ^(} r )-p(-y-; T ) -Plo 


v 3 (t) = 3( -p( };t)- p(-yh r I + 


(H 
2p (i ;T ) 
(H)' 


U(t) = 3 (-P ();r)+2p(^-tt;r J - p 

These two ways of writi ng the potentials make the action of dualities manifest. For instance, the transfor¬ 
mation properties (B.31 show that under S-duality, p(|, r) becomes p(|, z ^) = r 2 p(|, r) while p( r y^, r) 
becomes t 2 p{ t ^/t), so that V 2 and V 3 are S'-dual, et cetera. The result can also be written using perhaps 
more familiar modular forms 


V 2 {q) = - 6 ? T 2 E 2 , 2 {q) 

V 3 {q) = (2E 2 j 2 (g) - 3 ^( g )) 

V 4 (q) = (2E 2 j 2 (g) +30|(g)) . 

The action of T-duality is again clear from these expressions. For S'-duality it is slightly more intricate. 
Given that E 2 ^{q) = —0 2 (Q 2 ) ~ ^ 3 (? 2 )i it relies on the identities 

204(2r)+20 2 4 (2r) + 304( r ) = - 0 *( t / 2 ) + 20 4 (r/2) 

0 4 (r/2)+0 4 (r/2)-60|(r) = -40 4 (2 t) - 40 2 4 (2r) + 60 4 (r), 

for S-duality between extrema 2 and 3, and self-S-duality for extremum 4, respectively. 

11 There exist algorithms to find the generators. These are for instance implemented in Sage. 
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2.7.3 The quadruplet 

We move on to discuss the extremal values of the potential in the quadruplet. We can arrive at the 
following closed form for the potential in extremum 6 : 

V 6 (q) = -24n 2 (-^E 2t3 (q) + ( 77 (g ) 3 + Vrj(q 9 ) 3 )v{q 3 ) 2 /v{q) + 3(??(g 3 ) 3 A7(g)) 2 ) • 

Note that this can alternatively be written as 

V 6 (g) = -24t r 2 (g 0 (g) + g 1/3 gi(g) + 3 q 2/3 g 2 (q)), 

where the gi are functions that can be expanded into series with only integer powers of q (and the three 
summands in this expression correspond to the same summands in the expression above). Thus we know 
how the operation r —> t + 1 acts on the extremum, and it generates two other extrema, whose potential 
we also know exactly. These are extrema 7 and 8 : 

V 7 (q) = -247r 2 ( 5o (g) + e 2i7r/3 g 1/3 (/i(g) + 3e~ 2i7r/3 q 2/3 g 2 (q)) 

V»(q) = -247T 2 (g 0 (g) + e" 2i,r/3 g 1/3 gi(g) + 3e 2i7r,3 q 2/3 g 2 (q)) ■ 


The potential for the extremum 5 is: 


V 5 (q) = ~3TT 2 E 2 , 3 (q) ■ 

In the basis {5,6, 7, 8 } the matrices for S- and T-dualities are : 


f 1 

0 

0 



(° 

1 

0 

°\ 

0 

0 

0 

1 

s = 

1 

0 

0 

0 

0 

1 

0 

0 


0 

0 

0 

1 

V° 

0 

1 

0) 



0 

1 

0/ 


We can also apply the same method as above. The generators of T(3) are all trivial in this basis. Thus 
the potentials are weight 2 modular forms of this congruence subgroup. The latter form a 3-dimensional 
space, generated by the zero-sum linear combinations of the 4 Eisenstein series associated to the order 3 
vectors in (Z 3) 2 (there are 8 such vectors, but the Eisenstein series are invariant under v —> —v, leaving 
only 4 distinct functions, see appendix). We find 


V 5 = y[3G ; 

Ve = ?(-G 


* 2,3 


* 2,3 


V7 = — ( -G 2 ,3 


= 


27 


-G 


2,3 


— G 2> 3 
+ 3G 2i 3 
~ G 2 ,3 
~ G 2 ,3 


~ G 2> 3 
— G 2i 3 
+ 3G 2i 3 

— G 2 ,3 


G 2j 3 

— G 2j 3 

— G 2j 3 
‘ 3G 2j 3 


or alternatively, 


UM = 5(3p(j;r)-p(I ! T)-t>( 1 
um = 5(-p(y)+sp({{-p( 

(y T ) +3p 


V 7 (t) = 


-p ( 3 ;r )- P 


3 

T + 1 

3 

T+l 


;t -p 


t + 2 


-;r 


;t -p 


t + 2 


t 


;t - p 


v * (t) = ^(~ P (^ ;r ) _P ^ ;T )~ P (^ ;r ) +3p ( 


t + 2 
3 

t + 2 


t 


The dualities act on the vectors characterising the modular forms as follows 


o' 


o' 





1 


1 





' 1 ' 


Y 


'1 


' 1 ' 

0 

-+ 

1 

-+ 

2 

-+ 

0 


14 



s 


0 

1 

1 

1 




1 

0 

1 

2 


This reproduces the action of the dualities on the associated extrema. Thus, while the pattern of the 
positions of the extrema is non-linear, the arguments of the values of the potential at certain extrema do 
provide a linear realisation of the duality group. 

Finally, we note that triality generates three copies of the triplet as well as of the quadruplet. Indeed, 
each of these extrema is left invariant by a Z 2 subgroup of S 3 (as described in appendix D.l). 

Up to now, we have discussed the singlet, triplet and quadruplet whose duality diagrams are sum¬ 
marised in figure [5] 


2.7.4 The duodecuplet and a point of monodromy 

In the multiplet of size twelve, also depicted in figure [3j a new feature appears. We find that the extrema 
exhibit a monodromy around a point in the interior of the fundamental domain of the parameter r. Thus, 
to be able to describe the multiplet structure in this case we must first discuss the monodromy. 


The point of monodromy 

We find a single point in the interior of the fundamental domain around which there is monodromy 
amongst extrema. It is possible to determine this point numerical!}}* 2 ] and its value is close to tm ~ 
2.41558*. In particular, the extrema 13 and 16 are exchanged when we follow a loop in the r-plane that 
closely circles the value tm ■ Moreover, using the geometry of the positions of the extrema 13 and 16, one 
can show that tm is a solution of the system of equations 

f p(z;r) 2 + p(z — w 3 ;r ) 2 + p(2z — w 3 ;r ) 2 = ^S 4 (r) * 2 33 ) 

I 2 p'(~; T ) + 2 p '(2 - w 3 ; r) + p'(2z - w 3 ; r) = 0 , 

where UJ 3 = uq + w 2 , which gives the numerical result 

t m = 2.415576987549484510777262081474158860468152563579077460...*. 


Using the large accuracy of the value of the point of monodromy tm, we find the corresponding rational 
Klein invariant (with the normalisation (C. 1)): 


3{tm) 


488095744 

125 


= 1728 x 


7626496 

3375 


This can be considered as an exact statement - the uncertainty is as low as 10 200 . Elliptic curves with 
rational Klein invariant have interesting arithmetic properties (see e.g. [32]). 


The extended duality group 

We can add the monodromy group to the set of generators S and T that act on our vector of extrema. 
The resulting diagram of dualities then becomes the one in figure [3] The generators satisfy the relations: 

• S 2 = M 2 = 1 and T 6 = 1, while (TM ) 8 = 1 

• SM = MS 


• ( MST ) 3 = 1. 


Once we are underneath the point of monodromy in the canonical fundamental domain, the matrix A IT 
plays the role usually taken by the matrix T in SX(2,Z). In particular, relations like (ST ) 3 = 1 implied 
by the geometry of the fundamental domain of the modular group take on the form ( SMT ) 3 = 1, et 
cetera. Triality leaves each extremum invariant. 

we give terms in the Fourier expansion of the extremal values of the potential in 


In appendix D.l 


the duodecuplet. We note that a consistency and exhaustivity check on all multiplets is provided by the 
fact that the sum of all extrema in a given multiplet of SL( 2, Z) has to be a weight 2 modular form. The 
check works out: the sum equals zero in each multiplet separately, as it must. An analytic understanding 
of the duodecuplet extrema remains desirable. 


12 The most immediate manifestation of the monodromy phenomenon can be seen as a symmetry breaking in the equi¬ 
librium positions for extrema 13 and 16 when moving on the imaginary axis across the point of monodromy tm (which is 
purely imaginary). Below this critical value, as can be seen in the diagrams drawn at r = i (in appendix 1 1). I [ i. the two 
extrema are exchanged by the Z 2 action .Y, ea — X t , while above the critical value, they are both invariant with respect to 
this action. This makes it possible to determine 2.41557 < Imrjvf < 2.41558. 
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Figure 3: The diagram of the action of dualities on the D 4 = so( 8 ) extrema. In red we exhibit the action 
of S-duality, in green, T-duality, and in dotted blue, the monodromy. 


2.8 The Dual Cases B 3 = so(7) and C 3 = sp( 6 ) 

2.8.1 Exact multiplets 

For the twisted elliptic integrable models associated to the dual Lie algebra root systems so(7) and sp( 6 ), 
we present our results succinctly. We have found 17 isolated extrema for each, and they are Langlands 
dual. We have therefore 34 extrema in total. We identified two quadruplets of the full duality group 
for which we found analytic expressions for the potential at the extrema. The list of the corresponding 
extrema is given in appendix |D.2| We find the following duality properties and analytic values for the 
extrema of the potential. The extrema labelled {1,2} have extremal values for the so(7) potential equal 
to V\ (r) and V 2 (t). From the diagram of dualities (figure [4]), we read off that these extremal values 
are modular forms of To(4) with weight 2. Moreover, Langlands duality then implies that Viv(2 r) and 
V 2 v(2 t) are also of that ilk. The space A 42 (Fo( 4 )) of these weight 2 forms has the two generators 

-E2, 2(r) = dl( 2 r) + 0 4 (2t) = 1/2 (0 4 (r) + 0 4 (r)) 

-E 2 A {t) = 30 3 (2t) = 3/4(0 3 (t) + 0 2 (r )) 2 . 

In terms of the generators, the extrema are: 

Vi (r) = tt 2 {-E 2 , 2 (t)- 2 E 2 ,a(t)) 

V 2 (t) = 7T 2 {—7E 2 2 (t) + 2£ , 2 j 4 (t)) 

V±v (2r) = 7r 2 (+E 2 ^ 2 (t) + 0£ , 2 j 4 (t)) 

F 2 v(2t) = 7 t 2 (- 2 £ 2i2 (t) + IE 2 A (t)) . 


For the other quadruplet under the full duality group, we have a similar story, with the happy ending: 

V 3 ( 2 t) = 7 r 2 / 6 (- 15 E 2 , 2 (r)+ 7 E 2 i 4 (r)) 

V 4 { 2 t) = 7r 2 /6(+9E 2j2 (r)-lE 2j4 (r)) 

V3V (r) = 8 ^ 2 / 3 (- 3 E 2 , 2 (r) + lE 2i4 (r)) 

F 4 v(r) = 8^ 2 /3(0E 2i2 (r)-lE 2 , 4 (r)) . 
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The action of Langlands S 2 duality as well as T-duality can be found explicitly using these exact expres¬ 
sions, for instance by exploiting properties of 9 functions. As an example, we note that the action of 
T-duality is summarised in the equalities: 

T 2 ,2^ t +0 = -2T 2 , 2 (t) + T 2j4 (t) 

£ 2,4 ( T + 2 ^ = — 3T 2 2 (t) + 2E 2 ^(t) . 


Moreover, on the extrema, the Langlands duality S 2 acts as 


1 

2^2 


Vi ^97) = Viv(t) + 3tt 2 T 2 , 2 (t), 


and similar relations hold for the other 5 2 -dual couples, as predicted by the duality formula (2.8). 


2.8.2 The duodecuplet, the quattuordecuplet and the points of monodromy 

We further identified a duodecuplet and a quattuordecuplet under the duality group (for a total of 
(4 + 4 + 12 + 14)/2 = 17 extrema for B 3 = so(7)). Sufficient data to reproduce them is provided in 
appendix |D.2| These multiplets exhibit points of monodromy, and the full duality diagram is captured 
in figure [4] It should be understood that we only represent points of monodromy that are inequivalent 
(where two monodromies are taken to be equivalent when they are equal up to conjugation by other 
elements of the duality group). For instance S 2 M T S 2 is the monodromy around —l/(2r). 

We draw attention to a few features of the diagram. There are 5 extrema that form a quintuplet under 
T-duality (around r = zoo), labelled 5, 6, 7,8,9. When we also turn around the point of monodromy, 
the quintuplet enhances to a septuplet. This is reminiscent of a feature of the duality diagram for the 
duodecuplet of so(8). 

Finally, we performed an exhaustivity check on the extrema by summing the extremal values of the 
potential. We founcp^l 


E 

ie 4i 

= -87 t 2 E 2 i2 (t) 

E v i(r) 

i 64 2 

= 27r 2 T 2j2 (r) 

E y *( r ) 

ie 12 

= —207r 2 T 2i2 (r) 

E w 

ie 14 

= 19-7T 2 E 2:2 (t) , 


showing again that the sum of potentials over every multiplet is a modular form of To (4). 

This concludes our systematic case-by-case discussion of the low rank B, C , D isolated extrema of 
(twisted) elliptic Calogero-Moser models. We finish the section with a few further remarks on general 
features of the problem of identifying isolated extrema. 


2.9 Limiting Behaviour 

We wish to make a remark on the limiting behaviour of the integrable models near an extremum. As 
an example, consider extremum number 7 for B 2 = so( 5) which has its extremal positions equal to a 
real number plus r/2. We can take the limit of the potential as r -9 zoo while keeping the difference 
between the extremal positions and r/2 fixed. The limit of the integrable system is then a (trigonometric) 
Sutherland system of type D 2 = so(4), and indeed, the real part of the extremal positions agrees with 
those of the Sutherland system. This is but one example of the limiting behaviour of the models near 
the extrema. 


2.10 Partial Results for Other Lie Algebras 

In this subsection, we discuss very partial results for some higher rank Lie algebras. We think of the 
elliptic integrable model as a perturbation of the Sutherland model, with trigonometric potential. The 

’We evaluated the sum of the extrema numerically at two different values of r to identify the linear combination of E 2 ,2 
and £ 2,4 that equals the sum. We can then perform arbitrary many numerical checks at other values of r, and these work 
out. 
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Figure 4: The diagram of dualities for so( 7) and sp(6) extrema. In red, we show the action of Langlands 
^-duality on the extrema, in green, T-duality, and in dotted blue, monodromies, with the corresponding 
approximate values of the points of monodromy r. As discussed in the text, monodromies relating sp(6 ) 
extrema exist but are not represented here as they are equivalent to those already depicted. 


18 













Figure 5: The positions of the monodromies (red dots) inside the fundamental domain of r 0 (4) (shaded) 


Sutherland model has a ground state with all particles sprinkled on the real circle. We can perturb 
this traditional ground state by turning on the elliptic deformation by powers of the small parameter 9 , 
and follow the ground state under perturbation. In this way, we can reconstruct the extremum of the 
complexified elliptic potential associated to the Sutherland extremum on the real line. To take the limit 
from the elliptic integrable system towards the Sutherland model, it is sufficient to use the expansion 
formula: 


p(^;wi,w 2 ) 


2 2 




E 


nq n 
1 -q n 


mrx 

cos-, 

w 1 


(2.34) 


valid when the imaginary part of the modular parameter r is sufficiently large. The first term in the 
formula (2.34) is constant from the perspective of the integrable system dynamics, while the second 
term gives rise to the leading Sutherland potential. The minimum at the equilibrium of the Sutherland 
potential on the real line can be computed analytically na - it is related to the norm of the Weyl vector 
of the Lie algebra. The positions of the equilibria are given in terms of zeroes of the Jacobi polynomials. 
We can perform perturbation theory around these extrema (numerically), and we find the following series 
in q for the potential at perturbed Sutherland extrema, for various gauge algebras: 


^so( 5) 


7r 
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^so( 6) 
7T 2 

Vso( 7) 


V„, 


'( 8 ) 


K, 


'(9) 


26 112g 8656g 2 3921289 s 5011568 9 4 3191174729 s 

~3~ + ”3“ + 81 + 2187 + 19683 + 1594323 

122360348809 s 112088197760<? 7 

+ 43046721 + 387420489 + "' 

8 + 64 q+ 192 q 2 + 25 69 s + 192g 4 + 384g 5 + 768 q 6 + 512g 7 + 192q 8 + ... 

or 408g 153816 g 2 237305289 s 6103562136q 4 6633461285289 s 

+ + 625 + 78125 + 9765625 + 1220703125 

i 1293168139431369 s | 10819167546478272 9 7 
+ 152587890625 + 19073486328125 + 

0 576g 212544 9 2 395389449 s 66182636169 4 9098716296969 5 

+ ”5” + 625 + 78125 + 9765625 + 1220703125 

i 1714036086397449 s | 81126438184719369 7 | 10878191192254884489 s 
+ 152587890625 + 19073486328125 + 2384185791015625 

164 9929 51337289 2 23058446089 s 1689027994381129 4 

+ ~Y~ + 12005 + 4117715 + 176547030625 
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k 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

N 

3 3 

1 

5 3 

5 3 

7 3 5 2 

7 3 

3 a 7 2 

3 B 7 3 

3 1 7 1 11 3 

3 2 7 4 11 3 


Table 1: The integer N for gauge algebra so(fc) rendering the Fourier expansion integral 

11307570247017024g 5 640315787843154194816g 6 1106383118191321793331968g 7 

+ 12111126300875 + 370903242964296875 + 890538686357276796875 

69929754265259380435436903968g 8 17683503230173163609024329488224 9 9 

+ 38181846177568242666015625 + 13096373238905907234443359375 + 

Ko(io) 160 1280(7 1303808q 2 616518656g 3 365560247552g 4 

7r 2 ~ ~~T + 7 + 2401 + 823543 + 282475249 

101140172889600 9 5 9869502718168064<7 e 

+ 96889010407 + 4747561509943 

18401127697466238976(7 7 6582207315175560008960g® 

+ 11398895185373143 + 3909821048582988049 + 

Ko(ii) 305 1960q 30141880 9 2 29034410080g 3 4243088924219480 9 4 

7r 2 " ~~T + 9 + 45927 + 33480783 + 2790589782267 

7560807432828504560<? 5 4158609757083162994374880(7 6 

+ 6103019853817929 + 1526041805387611692663 

i 96348742286518866720674240g 7 | 304885265038041162579660724924120g® 

+ 52975451244169948759587 + 92724468600756742242419154123 + 

Ko(i 2 ) 800q 4055200g 2 1335804800g 3 63808646477600g 4 

7r 2 “ 1 + ~3~ + 5103 + 1240029 + 34451725707 

42945633858692800g 5 6332155765834649948800<7 6 

+ 25115308040403 + 2093335809859549647 + 

We see that at least for some extrema, it is fairly straightforward to generate interesting data on the value 
of the potential at these extrema at higher rank. We note a first pattern, valid at the order to which we 
have worked, in both the rank of the gauge group, and the power of the modular parameter q. Table [l] 
gives the conjectured smallest integer N such that for gauge algebra g, the potential V g {Nq)/n 2 has a 
Fourier expansion with only integer coefficients in the following sense: the expansion can be written as 
7 io(r + niq + n, 2 d 2 + ■■■) where the n* are integers, and the first term r is rational. 

As an example of this pattern, let us quote the formula: 

66679200tt 2 ^ o(12) ^ x q ) = ^792 +q + 745U3q2 + 252572301828g 3 + 108583732036588599 g 4 

+25066769592690393853446g 5 + 11087973934403204342320752348g 6 
+1966652180387341854168182867614728(7 7 + ... 

As a final remark, we note that our numerical searches in this and previous subsections are far from 
exhausting the capabilities of present day computers. 

3 The Massive Vacua of J\f = 1* gauge theories 

In this section, we first briefly review properties of the infrared physics of the Af = 1* supersymmet¬ 
ric gauge theory, and then show how the data we gathered on elliptic integrable systems in section [ 2 ] 
elucidates the physics of this gauge theory further. We obtain the M = 1* gauge theory from J\f = 4 
supersymmetric Yang-Mills theory with gauge group G by adding three masses m, for the three chiral 
M = 1 supermultiplets. We can then go to the Coulomb branch of the gauge theory, and compactify the 
theory on a circle mm. Two massless scalars remain in the theory for each U(\) in the unbroken U{l) r 
gauge group, namely the Wilson lines <j> = f sl A^dx^ and the scalar duals a of the photons. Since there 
are no fields in the theory which are charged under the center of the gauge group, we may choose the 
gauge group such that we allow for gauge transformations that twist around the circle by an element of 
the center. This lends a periodicity to the Wilson line under shifts taking values in the co-weight lattice 
P w . This reasoning corresponds to a choice of gauge group G = G/C where G is the universal cover, and 
C its center! 44 ] 

14 Note that the dual theory to the one with gauge group G/C has gauge group G, for a simply laced group. The two 
scalars are interchanged under S-duality. Thus, the duality symmetries mix various global choices of gauge groups. Duality 
also acts on the twist direction of the twisted elliptic potential. 
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The gauge theory compactified on a circle gets non-perturbative superpotential contributions from 
magnetic monopole configurations whose charges take values in the dual root lattice Q v . The scalar duals 
of the photons have as a result a smallest possible periodicity equal to the weight lattice P. We choose 
to classify extrema of the superpotential with respect to these identifications. We should mention that 
other choices will be physically relevant. Since in deriving the effective superpotential we compactified 
the theory on R 3 x S 1 , the resulting effective theory is influenced by the choice of the spectrum of line 
operators that probe the phases of our four-dimensional theory [3S£ :5J SI]. These determine the set of 
allowed monopole operators in three dimensions, and this set may be larger than the collection allowed by 
the minimal periodicity relation chosen above. Depending on the choice of the spectrum of line operators, 
this can lead to an increase of the number of inequivalent solutions, and therefore to an increase in the 
Witten index. This was analyzed carefully in [34| 33 • 

We have identified the shift symmetries acting on the Wilson line and the dual photon. We further 
divide out the configuration space by the Weyl group, which is the remnant of gauge invariance. This 
classification of supersymmetric vacua agrees with the classification we did in section [2] in the elliptic 
integrable systems, on the condition that we identify the 1 x 2 direction with the Wilson line. 

Our Af = 1* theory is a deformation of Af = 4 theory, and it inherits some of its properties. In 
particular, the electric-magnetic duality group of AA = 4 gauge theories in four dimensions [351CE1I3B! 
plays a crucial role. The duality symmetry was determined to be the group SL( 2,Z) for simply laced 
gauge groups and T 0 (4) for the B and C type gauge groups 39] 43s 31] ■ Moreover, the S 2 generator of 
the Hecke group exchanges the B and C type systems. An infrared counterpart to these duality groups 
are present in our integrable systems, which allow for a (generalized) duality group action on the infrared 
modular parameter r [26j , inherited after mass deformation from the AT = 4 duality. Note in particular 
that the requirement of the B type and C type exchange is implemented in our integrable system by the 
Langlands duality we discussed in subsection |2.2| This duality provides a further consistency check on 
the relative weight of the short and long root contributions, fixed in m through consistency with the 
superpotential of the pure AT = 1 super Yang-Mills theory. 

In [14] . following the reasonings in m m n na, an exact effective superpotential for AT = 1* was 
proposed for any gauge group, equal to the potential of the twisted elliptic Calogero-Moser model. The 
arguments were based on holomorphy, uniqueness of the deformation from J\f = 2*, the form of non- 
perturbative contributions, and integrability. We have added to these reasonings the test of S-duality in 
subsection 2.2 We wish to further strengthen the arguments for the superpotential by comparing the 
results for the exact quantum vacua for the theory on R 3 x S 1 with semi-classical results. 


3.1 Semi-classical Vacua 

A semi-classical analysis of the massive vacua of J\f = 1* on R 4 proceeds in several steps. First one solves 
the equations of motion for constant scalar field configurations which are equivalent to the statement that 
the three complex scalars satisfy a su( 2) algebra. The enumeration of inequivalent embeddings of su(2) in 
the gauge algebra then provides the set of classical solutions. In a second step, one analyzes the unbroken 
gauge group for each classical vacuum, and then counts the number of vacua that the corresponding pure 
J\f = 1 quantum theory gives rise to in the infrared (using e.g. the index calculation [42]). For gauge 
algebra su(n) the number of classical vacua was thus argued to be equal to the sum of the divisors of n 
[5], and this number coincides precisely with the number obtained from the exact superpotential [5] [TJ 
for the theory on R 3 x S 1 (where one classifies vacua in the manner described above). For other gauge 
algebras, the semi-classical counting of vacua was performed in [43] . For gauge algebra so(n) it was 
argued to be: 

Z semi-class ( x , y ) = X + x + x ^y + ^ + Q x * + gfi ( 6 + y) + x * ( 7 + 3 y) + ^(15 + 2 y) 

+x & {26 + j/ 2 ) + £ 9 (31 + 5 y) + ... (3-1) 

where the power of x is equal to n and the power of y is the number of massless I7(l)’s in a given massless 
branch of vacua. Although we will concentrate on massive vacua, let us remark that the semi-classical 
counting of massless vacua may well be futile in the full quantum theory, where there may be a single 
manifold of massless vacua of given rank I5F (albeit with different branches). Thus, we will only further 
consider the semi-classical formula © for y = 0. For low rank then, the formula gives the following 
number of massive vacua, semi-classically: 

so(5) : sc 6 so(6) : sc 7 so(7) : sc 15 so( 8) : sc 26. (3.2) 

We know the result for so(6) = szt(4) to be in agreement with the number of massive vacua of the exact 
superpotential for the theory on R 3 x5 1 . We moreover have that this counting of so(2r +1) vacua agrees 
with the semi-classical counting of vacua of sp(2r) [45], both on R 4 . In the following, we compare the 
predictions for the number of vacua some low rank gauge theories on R 4 to the results we obtained for 
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the massive vacua coded in the superpotential on R 3 x 5 1 . To make the comparison, we need to go into 
a little more detail of the semi-classical analysis. 


3.2 Low Rank Case Studies of Quantum Vacua 

In this subsection, we compare the analysis of integrable system extrema to the semi-classical analysis 
of massive vacua of Af = 1* gauge theory on R 4 case by case. We will moreover refine the counting at 
some stages by taking into account the transformation properties of the vacua under remaining global 
symmetries. This will also be the occasion to interpret the many duality properties that we found for the 
integrable systems in section [2] We also briefly comment on the monodromies. 

To wrap up a loose end first, let us note that the minimal mass Mj of a given vacuum i can be 
computed using the equation 

Ml = min [Spec(M^Mk)\ , (3.3) 

where A4k is the matrix of second derivatives of the potential in vacuum k: 

v k>% 3 dXidX.j 


This clarifies the logic behind our definition of isolated extrema of the integrable system (see equation 
(2.12) and the corresponding footnote). 


The case so(5) 

Semi-classically, we expect six vacua for the gauge theory on R 4 . Let’s recall in a little more detail 
how this counting arises. We allow for various five-dimensional representations of su( 2) as vacuum 
expectation values for the three complex scalars of Af = 1*. Even-dimensional representations must 
appear in even numbers. They need to take values in the gauge Lie algebra, and we classify them up to 
gauge equivalence. One then finds the following allowed representations [43) - we indicate the dimensions 
of the su(2 ) representations, the unbroken part of the gauge group, and then the number of massive 
vacua they give rise to in the infrared: 


5 

3 + 1 + 1 
2 + 2 + 1 
l+l+l+l+l 


1 : 1 
so(2): 0 
sp{2 ) : 2 
so(5) : 3. 


(3.4) 


For instance, the 2 + 2 + 1 dimensional representation breaks the gauge algebra down to sp( 2). Classically, 
this gives rise to a pure Af = 1 theory with sp(2 ) gauge algebra at low energies, which gives rise to two 
quantum vacua. Summing all the resulting numbers of semi-classical vacua, we find six massive vacua in 
total. 

When we compare this analysis to the exact quantum vacua that we found for the so( 5) gauge theory 
on R 3 x S' 1 , we remark that we have a neat correspondence. In particular, there is one vacuum, on the 
real axis, that we can identify in the exact quantum regime as the fully Higgsed vacuum (corresponding 
to the 5-dimensional irreducible representation of su( 2) in the list ( |3.4| ) ). Its S-dual we interpret as a 
confining vacuum, and it is a triplet under the T-transformation, agreeing neatly with the so(5) confining 
vacua (corresponding to the trivial representation of su( 2) in (3.4)). We moreover found a doublet under 
T-transformation, again in agreement with the two vacua corresponding to the sp( 2) classical vacuum. 
Thus, at this level, we find excellent agreement. We note that the analysis of section [2] demonstrates 
that the six vacua are in a single SL( 2,Z) sextuplet and that their transformation properties are in 
correspondence with the transformation properties of sublattices of the torus lattice. Their (generalized) 
S-duality and T-duality properties are now entirely known. 

The exact analysis has revealed a seventh vacuum on R 3 x S 1 . Its origin lies in the massless vacuum on 
R 4 . After compactication, we can render the massless vacuum massive. Indeed, we can turn on a Wilson 
line that commutes with the semi-classical configuration for the adjoint scalars, and that simultaneously 
breaks the abelian gauge group. The necessary Wilson line is precisely the one we found in the seventh 
quantum vacuum. We have thus found its senri-classical origin. 

One can wonder whether our identification used for the dual of the photon (mentioned in the intro¬ 
duction to section [3]), and therefore of the parameterization of the Coulomb branch moduli space reduced 
the number of physical vacua on R 3 x S' 1 . Indeed, identifying our model as the one corresponding to gauge 
group SO(5) + (in the nomenclature of [34], :35j), leads to a doubling of the triplet in the senri-classical 
analysis, while the other multiplets remain unchanged. For the 1 + 2 + 2 senri-classical split, this is the 
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case because the commutant is a SU{2) C SO( 5) gauge group (corresponding to a long root in 50(5)), 
and thus the pure AT = 1 gauge theory gives rise to only a doublet of vacua upon compactification. 

In the integrable system, this more careful analysis corresponds to the rule that solutions can only be 
identified under shifts by 2wi (and not u>i). A look at the so(5) extrema in the diagrams in subsection 


2.5 shows that this relaxed equivalence relation adds precisely three vacua, namely each of the confining 

Thus, 


vacua (labelled 2, 3 and 4) obtains a partner, as expected from the analysis of pure M = 1 
in this more careful treatment, we increase the number of vacua by three on both sides of the analysis. 

We have computed the masses of the vacua. They are all roughly of the same order, and much above 
the accuracy of our numerical approximations, thus guaranteeing that our vacua are indeed massive. 
Moreover, for a given massive vacuum, the values of the masses are all approximately within a factor of 
100 from each other. Interesting patterns in the (ratios) of masses (of various vacua) exist - it should be 
fruitful to study them systematically. 


The case so(8) 


In the case of the gauge algebra so(8), we find a further set of subtleties. First, let’s compare the quantum 
vacua on R 3 x 5 1 to the semi-classical analysis on R 4 . The semi-classical analysis yields [33] : 


7+1 

H : t. s 

5 + 3 

H : 1 , 

5+1+1+1 

so( 3) : 2 

4 + 4 

sp(2) : 2* 

3+3+1+1 

so( 2) x so(2) : 0 

3+2+2+1 

sp( 2): 2 S 

+1+1+1+1+1 

so(5) : 3 

2 + 2 + 2 + 2 

sp( 4) : 3* 

+2+1+1+1+1 

sp( 2) x so(4) : 6 

+1+1+1+1+1 

so(8) : 6 S 


(3.5) 


for a total of 26 massive vacua. The semi-classical analysis was done under the assumption that the Z 2 
outer automorphism of so[2N) is a gauge symmetry [43] , in contrast to our analysis in section [ 2 ] If we 
adopt this point of view, we are left with a single Z 2 global symmetry, and we have indicated in the 
counting above whether a set of vacua are a singlet (s) or are conjugate (*) under that remaining Z 2 . 

Using this analysis, and the T-duality transformation properties of the integrable system extrema, 
we can partially match the list of semi-classical and quantum vacua on R 4 and R 3 x 5 1 respectively. 
The 6 under the T-duality group makes for a match between extrema 10,11,12,13,14,15 and the su( 2) 
representation l + l + l + l + l + l + l + l. These correspond to the confining vacua. The doublets which 
are conjugate under the remaining global Z 2 match extrema 3 and 4 (as well as their Z 2 reflections) to 
the representations 4 + 4 and 5 + 1 + 1 + 1. The conjugate triplets match 6,7,8 (as well as their Z 2 
reflections) onto 3 + 1 + 1 + 1 + 1 + 1 and 2 + 2 + 2 + 2. The smaller representations of the T-duality 
group are harder to match. We can still identify the Higgs vacuum with the extremum number 9, which 
lies on the real axis and which we can therefore follow all the way to weak coupling. For other extrema, 
it is harder to follow the change of effective description of the gauge theory dynamics from the ultraviolet 
to the infrared. There is again a seeming mismatch of one in the total number of vacua. The origin is 
the same as before: one extra quantum vacuum arises from the massless vacuum in R 4 by turning on the 
appropriate Wilson line after compactification on S 1 . 

Of course, our modular analysis of extrema again obtains a gauge theory interpretation. Recall that 
we found a singlet, triplet and quadruplet under the modular group. The modular group plays the role 
of a generalized duality group |26], acting on the effective gauge coupling in the infrared. 

Note that we also found a more surprising feature: a new duality group, with a generator corresponding 
to a monodromy around a point in the fundamental domain of the effective coupling that we used to 
describe our theory. We found a duodecuplet of vacua transforming under this new duality group. It could 
be very interesting to understand this group better in terms of gauge theory physics, or, as associated to 
the choice of parameterisation in the infrared. 

Again, the masses all lie very amply above our numerical accuracy, such that we can claim that we 
indeed identified massive vacua. Masses are again within a factor of 100 or so from each other, and 
exhibit interesting patterns that could be explored. 
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The cases so(7) and sp{ 6) 

For gauge algebra so( 7), the semi-classical analysis on R 4 predicts fifteen massive vacua, that arise as 
follows [HI To): 







7 

1 : 1 




5 + 

1 

+ 

1 

so(2) 

: 0 



3 + 

3 

+ 

1 

so(2) 

: 0 



3 + 

2 

+ 

2 

sp{ 2) 

: 2 

3 + 1 

+ 

1 + 

1 

+ 

1 

so( 4) 

: 3 

2 + 2 

+ 

1 + 

1 

+ 

1 

sp( 2) 

x so(3) 

l+l+l+l 

+ 

1 + 

1 

+ 

1 

*o(7) 

: 5. 


In the quantum theory on R 3 x S 1 , we do find a quintuplet under T-duality associated to the confining 
vacuum on the imaginary axis, dual to the Higgs vacuum on the real axis, near weak effective coupling. It 
enhances to a septuplet under T-duality at stronger effective coupling. From our analysis of the quantum 
theory on R 3 x S 1 we see that the theory permits two more quantum vacua, labelled 1 and 2. Again, 
we checked explicitly that these arise from turning on Wilson lines in the massless vacua on R 4 . On the 
sp( 6) side of the duality, the two extra vacua l v and 2 V arise from an unbroken sp( 2) gauge group (after 
breaking the abelian group that commutes with the 2 + 2 + 1 + 1 representation). We note again that 
the multiplet structure under T-duality is blurred at strong effective coupling. 

3.3 Tensionless Domain Walls, Colliding Quantum Vacua and Masslessness 

The point in the fundamental domain around which we have found a monodromy in the case of the 
so(8) gauge algebra, corresponds to a point at which two massive vacua have equal superpotential. At 
this point, a supersymmetric domain wall between the vacua becomes tensionless |46| l47| . The physics 
associated to such a situation is hard to discuss in detail, because of the difficulty of controlling the 
Kahler potential in gauge theories with Af = 1 supersymmetry only. Explorations of the physics in this 
regime can be found in [481 HU 141)1 . We note in particular that in a mass and cubically deformed Af = 1 
U(N) theory in [481 an extension of the Zjv action associated to shifts in the 0 angle of the gauge 
theory to 'Lin was observed due to the presence of a point of monodromy in an effective coupling. The 
T-operation (shifting the 9 angle of the gauge theory) in our situation is also crucially influenced by the 
presence of the point of monodromy : above the point of monodromy (at weak effective coupling), we 
find a Zjv _2 action, while below (at strong effective coupling), we find a Zjy action (for the case N = 8 
as well as for the case of N = 7). We also note that the collision of the extrema of the superpotential 
indicates the existence of an effectively massless excitation since there will be a zero mode for the matrix 
of second derivatives. The physics, or at least the properties of the effective description, seem close to 
the discussion in e.g. [JH]. It would be interesting to elucidate this point further. 


4 Conclusions and Open Problems 

We studied the isolated extrema of complexified elliptic Calogero-Moser models, and encountered a 
plethora of beautiful phenomena. The values of the integrable interparticle potential at the extrema are 
true vector-valued modular forms in some cases, allowing for an analytic determination of the extrema in 
terms of modular forms of congruence subgroups of the modular group. This gives rise to webs of extrema 
that form representations under the duality group of the model. The latter can either be a modular or a 
Hecke group. A more intricate phenomenon is the appearance of monodromies amongst a second class of 
extrema as we loop around a point in the fundamental domain of the modular group. The duality group 
is then enlarged to include the monodromy generator. We determined the action of these generators on 
extrema. Moreover, we provided a wealth of Fourier coefficients of the extremal potential. These analyses 
can be viewed as a considerable widening of the observation of the integrality of observables in equilibria 
of integrable systems. 

Secondly, we interpreted the results on extrema of Calogero-Moser systems in terms of mass-deformed 
Af = 4 supersymmetric gauge theory in four dimensions. We compared our results based on a low-energy 
effective action for the quantum theory on R 3 x S 1 to semi-classical predictions for the theory. The total 
number of quantum vacua matched the number resulting from the semi-classical analysis, provided we 
took into account massive vacua that originate in massless vacua on R 4 . A Wilson line on the circle can 
commute with semi-classical expectation values for the adjoint scalars, yet break the remaining abelian 
factors in the gauge group to give rise to massive quantum vacua, of either Higgs or confining type. We 
note that the appearance of extra vacua after compactification is typical of the B, C and D series. The 
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compactified theory manifestly differs from the theory on M 4 . We also performed a more refined matching 
of quantum vacua in certain cases, thus providing further evidence that the superpotential gives a correct 
description of the physics of AT = 1* theory on M 3 x S 1 . Furthermore, we noted that the precise multiplet 
structures in the quantum theory showed surprising features, including monodromy properties of the 
quantum vacua. 

It should be clear that we only scratched the surface of this intriguing domain at the intersection of 
integrable systems, modularity and gauge theory. The new features of the extrema that we laid bare in 
the B , C, D-type integrable models (compared to the A -type theories) prompts the question of the generic 
counting of the extrema, the relevant duality group and modular structure (including monodromies) as 
well as their representation and number theoretic content. Clearly, our analysis begs to be extended to 
exceptional algebras of low rank, to higher rank root systems, to models with different choices of coupling 
constants, as well as to the integrable generalizations of the elliptic Calogero-Moser models, including for 
instance the Ruijsenaars model with spin. Moreover, our study can be extended to other observables, 
like the ratio of the frequencies of fluctuations. It would also be interesting to attempt to characterize 
the positions of the extrema through e.g. a generalization of zeroes of orthogonal polynomials 125;. All 
indications are that similarly intriguing phenomena as the ones we uncovered will appear in this broader 
field. 

In gauge theory, one would like to analyze more closely the Seiberg-Witten curves of the Af — 2* 
theories that underlie our models, and in particular, locate the points in Coulomb moduli space where 
the curve develops a number of nodes equal to the rank of the gauge group, and where the vanishing 
cycles are mutually local (indicating the existence of massive vacua after mass deformation). The Seiberg- 
Witten curves are defined by equations of higher order, rendering this analysis harder than in the cases 
treated in detail so far [5j. 

One would also like to have access to the large rank generalization of our results, to connect to 
holographic dual backgrounds with orientifold planes mmm- For these purposes it might suffice to 
have access to the large rank generalization of a Higgs and confining vacuum, which one may hope to 
characterize analytically. It would also be useful to perform a more careful analysis of the discrete choices 
of gauge groups and line operators in our model [34' L35J, and to classify various supersymmetric vacua 
further [33) (M], e -g- by understanding a single phase, then chasing it through the duality chains. 

Finally, we already noted in passing that the branches of massless vacua predicted by the semi-classical 
analysis may turn out to be connected in the quantum theory, giving rise to a single massless vacuum 
manifold, consisting of branches that can be characterized by differing algebraic equations [511 . Our 
numerical explorations up to now are consistent with the fact that all massless vacua have the same value 
of the superpotential. It would be interesting to clarify the structure of these vacuum manifolds further 
for the models at hand. 
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A Lie Algebra 

We briefly review Lie algebra concepts that are useful to us in discussing the symmetries of both the 
integrable systems and gauge theories we discuss in the bulk of the paper. See e.g. [24] for a detailed 
exposition of the following facts. Let us consider a (compact simple) Lie group G with maximal torus 
T. They have corresponding tangent algebras g and t. We can then identify T as a linear group, and its 
space of characters x(T) is i n bijection with a lattice in the space t*(R) dual to the tangent algebra t, and 
defined over the real numbers R. To the Lie algebra, we can associate its space of weights in the adjoint 
representation, which is the set of roots A. Again, these roots are elements of the Euclidean space t*(R). 
The space t(R) comes equipped with a non-degenerate scalar product, which we will denote (•,•)• This 
scalar product allows us to identify a function A on the space t(R) with an element u\ of the space t(R) 
through the relation: 


A(z) = (u x ,x), (A.l) 

valid for all elements x of t(R). We will occasionally abuse notation and write A(x) = (X,x), and 
also = (A, A'), which defines a dual scalar product on t*(R). The bijection between the space 

generated by the roots and its dual allows us to define the dual roots (i.e. the co-roots) through the 
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relation: 


2 u a 
(a, a) ' 


a 


v 


(A.2) 


The root lattice Q is the lattice generated by the roots. Any set of simple roots a* spans the space t*. 
The weight lattice P also sits inside t* and is defined to be generated by a basis 7T ; - such that: 


r, {ctj, TTj) 


(A.3) 


for all z and j that run from 1 to the rank of the group G. We moreover define the dual root lattice Q v 
to be the lattice generated by the dual roots, and the dual weight lattice P v to be the weight lattice 
corresponding to the dual root lattice. The dual of the lattice generated by the characters of a given 
group G will be denoted t(Z). We have the following properties. The center C(G) of the group G is given 
by: 


C(G) = P^/t(Z) = X (T)/Q. (A.4) 

Moreover, when G is simply connected it is equal to its universal cover G. We then have that the space of 
characters is bijective to the whole of the weight lattice x(T) = P, and that t(Z) = Q v , such that C(G) 
is maximal and C{G) = P/Q = P v /Q w . The group with minimal center C(G) = 1 is the universal cover 
G divided by its center C{G). In this case we have that the set of weights is the set of roots x(T) = Q 
and that t(Z) = P v . 

Our definitions imply that the fundamental weights it)- that generate the dual weight lattice P v 
satisfy: 


(ttY, u ai ) = Si- 


(A.5) 


and therefore that: 


a(X) = (u a ,X) = (a, X) 


(A.6) 


is integer for X in the dual weight lattice, i.e. for X a co-weight. 

We summarize inclusions and dualities in the diagram below. The arrows indicate that the lattices 
are dual, i.e. that the contractions give integers. 


D 


p 

3 X(T) D 


X 

g v = p* 

c t(Z) c 


Q 

l 

P v = Q* 


c 


Weyl group, the outer automorphisms, the dual Coxeter number and the center of the universal covering 
group corresponding to the classical Lie algebras: 


Algebra 

Weyl group 

Outer Automorphisms 

Dual Coxeter number 

Center Univ. Cover 

A r , r > 1 

Sr-\- 1 

Z2 

r + 1 

Z r +i 


Z2 

1 

2 

Z2 

B r 

S r K Z^ 

1 

2 r- 1 

Z2 

C r 

S r K Z^ 

1 

r +1 

Z2 

D ri odd r 

S r K Z T ~ X 

Z2 

2i — 2 

7 L/± 

D r , even r > 4 

S r K Z^" 1 

Z2 

2j— 2 

^2 ^ ^2 

P 4 

S4 K Z| 

£3 

6 

^2 X ^2 


Table 2: Lie Algebra Data 


B Elliptic Functions 

Our conventions for the elliptic Weierstrass function are: 

p(s;wi,W2) = J2 

(m,n)^( 0,0) 

p( z ;r) = Y, 


1 


(m,n)^( 0,0) 


(x + 2muji + 2rvjj2) 2 (2mwi + 2nuj2) 2 

1 1 


(z + m + nr) 2 (to + nr) 5 


(B.l) 
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which entails the equality 


p(z-,T = u 2 /uj i ) = 4ujf p(2cv 1 z;cv 1 ,uj 2 ) ■ (B.2) 

We impose the convention that Q(ui 2 /uji) = 3(r) > 0. The Weierstrass function is a Jacobi form of level 
2 and index 0 : 


P 


ar + b 


ct + d ct + d 


(ct + d) 2 p{z;r). 


(B.3) 


It has the following expansion for large imaginary part of r: 


p(x;uj i,w 2 ) 


jE 2 (q) + —-o esc 2 


12w 2 “ 4 uif 




E 


nq n 
1 -q n 


nnx 

cos-. 

Wl 


(B.4) 


For the twisted Weierstrass functions, we can derive the equalities: 

UJ 7T^ id <jJ 

p(x;wi,w 2 ) + p(x +wi;wi,w 2 ) = p(x; —, w 2 ) + —^(2F1 2 (2 — )-.E 2 (—)) 

p(x; wi, w 2 ) + p(x + w 2 ; wi, w 2 ) = p(x; wi, - ^(Fl 2 ( —) - E 2 (^-)). (B.5) 

Z Z £(jJ\ 

These can be proven using the definition of the Weierstrass function p, as well as the definition of the 
second Eisenstein series E 2 . 


C Modular Forms 

We present a compendium of modular forms that we put to use in the bulk of our paper. 

C.l Theta and Eta Functions 

We first fix our conventions for the theta-functions with characteristics: 


8 


(r) = exp [in(n + a) 2 r + 2nil3(n + a)] . 


n£Z 


Particular examples of these theta-functions include: 

"l" 


2 ( 9 ) = e g 


(q) = 2 q 1/s + 2 q 9/s + 2 q 25/8 + 2 g 49/8 + ... 

03(q) = 0 Jj (q) = l + 2q 1 / 2 + 2q 2 + 2q Q ' 2 + 2q* + ... 

8i(q) = 8 ? (q) = 1 — 2q 1 / 2 + 2q 2 — 2q 9 ^ 2 + 2q 8 + ... 
. 2 . 

We also make use of the Dedekind eta-function: 

00 

ri(q) = q 1 / 2 i \{(l-q n ), 

n—1 

and the Klein invariant 


j(q) = 1728 


El(q) 


El(q) - El(q) q 


= —b 744 + 196884g + ... 


(C.l) 


C.2 Modular Forms and Sublattices 


In this subsection we recall how to find a basis of the space of modular forms of weight k for the congruence 
subgroup ((32]) : 


r w -{(: :)=() ;)- d 4 
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First we note that the cusps of T(N) can be identified with the pairs ±.v £ (Z/iVZ) 2 of order N. This 
makes it possible to count the number of such cusps : 


eoo(r(iV)) 



N = 2 
N> 3. 


For any congruence subgroup T the space of modular forms yVl fc (r) of weight k decomposes into the 
subspace of cusp forms and the Eisenstein space: M k {X) = <Sfc(r) ® ffc(T). For N = 2 we have 


dim5 2 (r(2)) = 0 


and for N > 3, 


dim 6 > 2 (r(./V)) = 1 + 


N 2 (N — 6) 
24 




In particular, dimiS 2 (r(3)) = 0 and dim«S 2 (r(6)) = 1. 

We want an explicit basis of the Eisenstein space. For any vector v 
and for k > 3, we define the (non-normalized) Eisenstein series 


c 

d 


£ (Z/7VZ) 2 of order N, 


Gk,N M (t) = Gfc ,N 


M= E' 7^ 


v'=v(N ) 


( c't + d') k 1 


and for weight two 


G2,n [w] (t) = G 2 ,at 


c 

, 1 

(CT + d \ _ . . 

d 

{T) = w 

[ P { N ’ t )+ G2 ^\ 


where the primed sum runs over those non-vanishing vectors v' 



that equal v modulo N. One can 


show that the Fourier expansion of these functions in terms of q = e 2lnT is: 


G k . 


N 


(<?) = 


(—2 ni) 


’ V" ^ 


(n)q 


,n/N 


where 


Cfc-l.AT 


N 2 {k - 1)! ^ 

\ (n) = ^ sgn(m)m fe_1 exp (2m^\ 


m|n and ^ =c(N ) 


and 


Ca r(k) 


= E' 

d'=d (N) 


1 


This Fourier expansion is valid for all k > 2, including k = 2 which is the case we are mostly interested 
in. 

For k > 3, any set {Gfc^u]} with one v corresponding to each cusp of r(iV) represents a basis of 
the space £ k (T(N)) of Eisenstein series of weight k on r(AT) (and in particular dim£*.(r(./V)) = Coo). For 
the case k = 2 these statements have to be modified, because of the lack of modularity of the (ordinary) 
weight 2 Eisenstein series. It turns out that dim£ 2 (r(7V)) = — 1, and that £ 2 (T(N)) is the set of 

linear combinations of the {Gfc^M} (where v £ (Z/iVZ) 2 is of order N) whose coefficients sum to 00 

The Eisenstein series G k ,N M have good transformation properties under SL( 2, Z) for k > 3 and 
N £ {2,3} provided the vector v is transformed accordingly P’1 


1 

(cr + d) k 


Gk,N [«] 


/ ar + b\ 
\ct + d ) 




(t)- 


For k = 2, we have to take into account a non-holomorphic term, except for linear combinations where 
the sum of the coefficients vanishes, as is the case for the potentials considered in the bulk of the paper. 

15 Theorem 4.6.1 in m 

16 For generic N the relation between the normalized Eisenstein series, which enjoy these good transformation properties, 
and the series G & n [i>] is not simply a proportionality relation (see formula (4.5) in m ), but it is a simple rescaling for 
N = 2 and TV = 3. 
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Finally, we also define 


E 2 ,n(t) = E 2 (t)-NE 2 (Nt). (C.2) 

These are weight 2 modular forms of Tq(N). We use extensively the fact that Ad 2 (To (4)) has dimension 
2 and is generated by 

-E 2>2 {q) = 1 + 2 Aq + 24q 2 + 96q 3 + 2 Aq 4 + 144g 5 + ... 

-E 2 A (q) = 3 + 24q + 72q 2 + 96q 3 + 72q 4 + 144g 5 + ... 

We note the transformation property of the form E 2j2 under r —> —1/(2r): 

£ 2 , 2 (-1/(2t)) = E 2 (—l/ (2t)) — 2£’ 2 (—1/r) = —2t 2 E 22 (t) . (C.3) 


D The List of Extrema 

In this appendix, we list the extrema of the complexified (twisted) elliptic Calogero-Moser models with 
root systems D 4 = so( 8 ) and B 3 = so(7). We provide a few more details on how we obtained them, how 
to relate them through dualities, as well as Fourier expansions of the extremal potentials. 


D.l The List of Extrema for so(8) 

The strategy we used to find extrema boils down to finding all the minima (which are also zeros) of 
the (auxiliary, gauge theory) potential (2.12) with non vanishing mass (3.3) using a simple gradient 
algorithm with random initial conditions. Then we identify those configurations that are related by one 
of the symmetries we quotient by. This procedure is executed at a given value of r. Once the complete 
list of extrema is known, we can follow a given extremum along any curve in the r upper half plane, by 
adiabatically varying r. The T-dual extrema and the monodromies are obtained in this way, while the 
action of S'-duality is known exactly. We thus unfold the whole web of dualities. 

In order to determine the potential at the extrema, we first make use of our knowledge of T-duality, 
which dictates the Fourier expansion variable q 1,/n , where n is the smallest positive integer such that T n 
acts trivially on the extremum under consideration. Then we evaluate the extremal potential at many 
different values of r and find recursively the rational Fourier coefficients. 


D.1.1 The diagrams of the extrema 

In the diagrams that follow, the black dots represent the values of the components X,, i = 1,2, 3,4 at 
the extrema. The dark grey dots are images under the symmetries discussed in section |2.5[ 

In some of the diagrams, there are five black dots instead of four, reflecting the fact that they represent 
three extrema related by the global S3 symmetry. For every such extremum, one subgroup Z 2 C S3 acts 
trivially. One of the three extrema is obtained by choosing one of the circled black dots and the three 
ordinary black dots, a second one is obtained by choosing the other circled black dot and the three black 
dots, while the third is determined by the four small black dots. (The pale grey dots show the possible 
translations of this extremum by half-periods.) 

We note in passing that some exact information on the positioning of the extrema is available. For 
instance, for extremum number 9, some exact information on the positions is the following. At r —► ioo, 
the system reduces to the Sutherland system (with trigonometric potential). According to [17], the 
positions at equilibrium are related to the roots of a Jacobi polynomial. Explicitly in the case of D 4 , 
the polynomial is P 2 ^\y) = ^{y — l ) 2 + — 1) + 3, from which we deduce the positions X\ = 0, 

X 2j3 = arccos(±l/v / 5) and X4 = /. For r —> 0, the positions converge on X-\ = 0, X 2 = 1/6, X 2 = 1/3 
and X4 = J. This numerical convergence is slow. 

S-duality guarantees that the situation is similar for the extremum on the imaginary axis, with the 
two limits exchanged. Moreover, T-duality then acts in the r —► ioo limit as Xq —> Xq, X- t —>• X\ + 1/6, 
X 2 X 2 + 1/3, X 3 —> X 3 + 1/2. (These transformations are exact within the precision of the numerics.) 
This generates the 6 -cycle. Et cetera. 
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Extrema at r = i for so(8) 




• • 



Extremum 2 Extremum 3 Extremum 4 



• • 


Extremum 8 
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Extrema at r = i 




• • 

—i-#- 


Extremum 9 


Extremum 10 


Extremum 11 


• • 

• • 

• • 

• • 


• • 

• • 

• • 

—•-*- 


• • 

• • 

• • 

• • 


Extremum 12 


Extremum 13 


Extremum 14 


• • 

• • 

• • 

—#-#- 


• • 
• • 


• • • • 


Extremum 15 


Extremum 16 


Extremum 17 


• • 
• • 


• • 

• • 
• • 


• • 

-§-#- 


• • 

• • 

• • 

• • 

• • 

• • 

-§-t- 


• • 
• • 


• • 

-f-© 


Extremum 18 


Extremum 19 


Extremum 20 


D.1.2 The series for the so(8 ) extremal potentials 

We have been able to determine the g-expansions of the potentials in each extremum with great accuracy, 
in terms of functions with integer coefficients. For extrema 1 to 8, we gave the exact expression in section 
|2.7[ To list the series for the remaining extrema, we introduce 11 functions, for which we only reproduce 
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2 


1 - 


0 5 10 15 20 

Figure 6: The dots show the successive ratios of the coefficients of /n, and a line has been drawn, for 
comparison, at the value 1 /qm = e~ 2mTM . 


the first few coefficients - more can be obtained - : 

/i {q) = ^ - 467g + 45379<? 2 - 23993958092 9 3 - 44044347374301g 4 - 711960536580667762g 5 + ... 

f 2 (q) = 1 - 15172(7 + 51582918g 2 - 397077052296 9 3 + 5101142359347277g 4 + 94300056917523369780g 5 + ... 
f 3 (q) = gig +g + 369 9 2 + 68644 9 3 + 11490041g 4 + 1579638246g 5 + ... 

U(q) = 1 + 3096g + 1818378264g 2 + 2446348866170976g 3 + 4535490919062930456600g 4 + ... 
f 5 (q) = 1 - 142284g - 2825331513294g 2 - 110241726267588876840g 3 + ... 

f 6 (q) = 2 + 780960g + 18367562372664g 2 + 762875530342634406144 9 3 + ... 

f 7 (q) = 1 - 4478868(7 - 121113750523626g 2 - 5314750232983801186536g 3 
f s (q) = ±(14 + 79929712 q + 2425403175787968 q 2 + 111756708524847535116096(7 3 + ...) 
f 9 (q) = f (-37 - 489421748g - 16364614670173794g 2 - 787663906596039662206584g 3 + ...) 

f 10 (q) = 1 - 12264 q - 7273512936 9 2 - 9785395464683424g 3 - 18141963676251721826280g 4 + ... 

/u(V) = 1 + 110596g + 110757888006q 2 + 180011523750912008g 3 + 367762906594569664954381g 4 + ... 

The potentials then read 


Eg = 

144007t 2 / 3 (£) 

Vio+fc = 

-4tt 2 (16g) J '/ 6 exp (27 f 4+j (§) 

Vl6 = 

-37T 2 (/i 0 (g) - 72 y/qf n (q)) 

Vn = 

-37T 2 (/i 0 (g) + 72 y /qf 11 (q)) 

Vfg = 

-24^ 2 (75/!(g/15 3 ) + i^/3f 2 (q/15 3 )) 

v 20 = 

-24^(75/,(<?/15 3 ) - i^/3f 2 (q/ 15 3 )), 


where k = 0,..., 5. The last series Vis can then be deduced from the fact that the sum of all potentials 
in the duodecuplet vanishes. Note that the coefficients grow rapidly, preventing the functions above to 
be modular forms. The monodromy is responsible for this phenomenon, as can be confirmed by the 
estimation of the convergence radius given by the successive ratios of the coefficients (see figure [ 6 ]). 


D.2 The List of Extrema for so(7) and sp( 6) 

Finally, in the case of the algebras B 3 = so( 7) and C 3 = sp( 6 ), we only present diagrams of the extremal 
positions for the so(7) root system, since the corresponding extrema for sp( 6 ) can be found by Langlands 
duality. We use the same conventions as for the so( 5) figures. Additional data, like the data we presented 
for so( 8 ) in the previous section, can be found. 


32 








Extrema at r = i for so( 7) 



Extremum 4 

Extremum 5 

Extremum 6 

• 

• 

• • 

7 • • 

• 

• 

• 

O • • 8 

• 

• 

• 

• • 

• 

• • 

Extremum 7 

Extremum 8 

Extremum 9 



Extremum 10 


Extremum 11 


Extremum 12 





Extrema at r = i for so(7) 



Extremum 16 Extremum 17 
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